Reversible modulation and ultrafast dynamics of THz resonances 
in strongly photoexcited metamaterials 
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We demonstrate an ultrafast reversible modulation of resonant terahertz (THz) response in 
strongly photoexcited metamaterials. The transient spectral-temporal response of the dipole tran- 
sition ^1.6 THz exhibits a distinct non-monotonic variation as a function of pump fluence. The 
transition energy shift, strength, spectral width and density-dependent ultrafast relaxation manifest 
a remarkable re-emergence of the resonances after initial quenching. Our simulation, incorporating 
the first-order diffraction from the photoinduced transient grating, reproduces the salient features, 
providing a new avenue for designing nonlinear and frequency- agile THz modulators. 

PACS numbers: 78.67.Pt, 78.20.-e, 78.47.-p, 42.25.Bs 



Artificially subwavelength structured materials, so- 
called metamaterials [l|-l4(, attract strong current inter- 
est due to their exceptional properties, such as simul- 
taneously negative permittivity e(uj) and permeability 
/i(cj), which implies a negative refractive index [l45j 
not available in natural materials. In addition, artificial 
magnetism @, , super focusing (g, ^ and specifically- 
tailored structures responding to ultra-broadband elec- 
tromagnetic radiation, from gigahertz (GHz), terahertz 
(THz) to visible range 0, [Io|, set them among the most 
promising candidates for next generation optoelectronic 
devices and large-scale photonic functional systems. 

Recently, there has been significant interest in under- 
standing ultrafast THz responses of nonlinear metama- 
terials, e.g., split ring resonators (SRRs) patterned on 
substrates exhibiting substantial nonlinearities [ll|. The 
resulting components and schemes allow for construct- 
ing controllable active THz devices with dynamical tun- 
ability of amplitude and phase, which fills the demand- 
ing THz technology gap. Prior experiments using time- 
domain THz spectroscopy have demonstrated dynamic 
tuning of both magnetic 12l-[l7| and electric [HI dipole 
resonances of SRRs on semiconductors or superconduc- 
tors. More complex structures of coupled resonators on 
GaAs substrates have been proposed and experimentally 
studied, revealing frequency- agile THz modulation, e.g., 
a resonance shift with increasing photoexcitation [14|, [l6| . 
However, thus far all ultrafast studies of THz materials 
have been concentrated on the relatively low photoexci- 
tation, where a reduction or Quenching of the resonant 
absorptions is observed [l3|, [3 [l9| . Two outstanding is- 
sues are still poorly understood: first, there lack insights 
on how high density photoexcited carriers modify the 
dynamic responses of metamaterials; second, the decay 
pathway of photoexcited metamaterials after femtosec- 
ond (fs) excitation is not explicitly studied, and, par- 
ticularly, the temporal evolution of the transient states 



have only been understood as a density- independent re- 
laxation back to the equilibrium. 

In this letter, we demonstrate an ultrafast reversible 
modulation of resonant THz absorption in strongly pho- 
toexcited metamaterials by optical pump and THz probe 
spectroscopy. This process directly manifests itself via a 
remarkable reemergence of the originally quenched THz 
resonance above a crossover density N c ~4.7 x 10 16 cm -3 
after femtosecond photo-excitation. Increasing the ex- 
citation from below to above the transition density, we 
identify two distinctly different relaxation pathways of 
the THz resonance with opposite dynamics. A model cal- 
culation of transient spectra incorporating the first order 
diffraction mode from the photoinduced transient grating 
reproduces the salient features, which is further corrob- 
orated by their dependence on the unit cell lattice con- 
stant. The revealed scheme represents a relatively simple 
and generic approach to achieve nonlinear and frequency- 
agile functions in THz device without particularly com- 
plex structures. 

The samples used are double SRRs on GaAs, illus- 
trated in the inset of Fig. 1 (a), which are patterned 
6 jam copper rings on high resistivity GaAs of 630 jam 
thick. We mainly focus on two SRR samples with unit 
cell lattice constant of 50 jam and 45 /im, and the outer 
dimension of an individual SRR is 36 jam. Our optical 
pump and THz probe spectroscopy setup is driven by a 
1 kHz Tiisapphire regenerative amplifier with 40 fs pulse 
duration at 800 nm. One part of the output is used to 
excite the sample, while a small fraction is used to gen- 
erate and detect THz pulses [20| . Phase- locked THz field 
transients are used as a probe, which is generated via 
optical rectification and detected by electro-optic sam- 
pling in 1 mm ZnTe crystals (supplementary material). 
Copper apertures of 2.5 mm are placed in front of the 
sample and a GaAs reference to ensure a uniform illumi- 
nation and faithful comparison. The THz field transmis- 
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FIG. 1. (color online) (a): Transient transmission spectra 
T(uj) taken under two excitation nuences, 19 /xJ/cm 2 (cyan 
dashed) and 57 /iJ/cm 2 (red solid), respectively. Static trans- 
mission spectra without pump is plot together (black dotted) . 
The inset: the SRR unit cell of our sample and the polar- 
ization of the normally incident THz radiation. a=50 /im, 
b=36 /im, d=3 /xm, g=2 /xm and w=6 /im. (b): The non- 
monotonic pump fluence dependence (indicated by the gray 
arrow) of the resonant absorption ^1.63 THz for low (dashed 
lines) and high (solid lines) photoexcitation. Spectra In both 
panels are taken at fixed time delay of 20 ps. 



mains quenched. This salient feature is further corrob- 
orated in the detailed pump fluence dependence shown 
in Fig. 1(b). In the low excitation regime from 1.2 to 
19 /iJ/cm 2 , there is a clear weakening in the transition 
strength and broadening in the linewidth for the cui res- 
onance, with an almost complete quenching at pump flu- 
ence of 19 /iJ/cm 2 . Here the transmission is ~45% higher 
and the peak energy has a blue shift of ~60 GHz. In 
strong contrast to this, futher increasing photoexcitation 
from 56 /iJ/cm 2 to 112 /iJ/cm 2 results in a progressively 
pronounced resonance. The spectra clearly show a red 
shift of the resonance energy, ~20 GHz for the highest 
pump fluence used, a partial recovery of the transition 
strength and linewidth narrowing. This remarkable non- 
monotonic behavior and re-emerging resonance near the 
uji manifest a reversible and frequency- agile modulation 
of the resonant THz absorption by photoexcitation. This 
behavior has been seen experimentally for the first time, 
and we can explain below that for very high pump flu- 
ence the first-order diffraction will surpass the electric 
resonance behavior in the THz responses of photoexcited 
metamaterials. Fig. 1(b), the dips in the transmission 
above 56 /iJ/cm 2 are due to the first-order diffraction 
peaks, which will be discussed further below. 



sion coefficient T(uj) is obtained by the ratio between two 
Fourier transformed spectra from the transmitted THz 
probe pulses through the SRRs and bare GaAs. The 
transient THz signals are recorded at various time delays 
with respect to the pump pulse to obtain time-dependent 
T(uj). The whole setup is purged with dry-N2 gas. 

The static transmission spectra are shown in Fig. 1(a) 
(black dotted) with a THz probe normally incident on 
the sample plane. The E field of the probe is parallel to 
the gaps of SRRs (inset). There are two main THz reso- 
nances detected in this configuration at uj\ ~1.6 THz and 
uo2~ 0.5 THz, which originate from resonant electric pro- 
cesses consistent with the prior measurements [l3|, |l8| . 
The uo\ can be understood as the electric dipole reso- 
nance of the metallic bars of the SRRs, and the U02 is 
the magnetic dipole resonance induced by the circulating 
electric currents generated by the incident E field. 

The two THz resonances exhibit distinctly different 
pump fluence dependence, as shown in the transient THz 
transmission spectra at At =20 ps (Fig. 1(a)) for two 
pump fluences: 19 /iJ/cm 2 (cyan dashed) and 56 /iJ/cm 2 
(red solid). The 20 ps time-delay between the excitation 
and probe pulse is introduced to ensure a measure of 
quasi-steady state considering the transient carrier life- 
time in GaAs on the order of 1 ns. The 19 /iJ/cm 2 trace 
exhibits largely quenched resonances for both the uo\ and 
U02 after the photoexcitation. However, the 56 /iJ/cm 2 
trace clearly shows substantial recovery of the resonant 
THz absorption around cji, while the UO2 resonance re- 




FIG. 2. (color online) Time evolution of negative transmission 
spectra, -T(cj), in the photoexcited metamaterial of 50/xm 
lattice constant for (a) 4 /xJ/cm 2 and (b) 70 /xJ/cm 2 , respec- 
tively. These clearly show the opposite relaxation pathways 
for the weak and strong excitation regimes (see text). 



The weak and strong photoexcitations lead to signifi- 
cantly different relaxation pathways back to the equilib- 
rium, as shown in the time evolution of T(cj) spectra fol- 
lowing photoexcitation in Figs. 2(a) and 2(b). For pump 
fluence at 4 /iJ/cm 2 , the THz absorption ~uj\ resonance 
exhibit an initial photoinduced reduction, immediately 
followed by a relaxation back to the original vaule on a 
ns time scale. Interestingly enough, opposite to those at 
the low excitation, the 70 /iJ/cm 2 dynamics clearly shows 
a further reduction of the THz resonence ~1.6 THz fol- 
lowing the initial decrease. We emphasize two key as- 
pects of this observation: (i) at the large time delay of 
Ar =800 ps, the THz resonance ~1.6 THz is almost com- 
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FIG. 3. (color online) (a) The extracted transient carrier con- 
centration as function of time delay. The left and right gray 
areas illustrate a crossover from the high to low excitation 
regime with a crossover density N c (black square). The inset 
plots the extracted transient carrier concentration as function 
of pump fluence at Ar =20 ps, marked together with N c cor- 
responding to the density excited by 19 /xJ/cm 2 . (b)-(e): The 
transient transmission spectra T(cj) at various time delays 
following 70 /iJ/cm 2 photoexcitation. 



pletely quenched for the 70 /iJ/cm 2 excitation, although 
a clear recovery of the resonance has been seen for the 
4 /iJ/cm 2 case; (ii) the pump-fluence-dependent relax- 
ation dynamics of the resonance indicates another way 
to achieve nonlinear and reversible THz modulation by 
tuning the time delay. 

Ultrafast photoexcitation strongly alters the electronic 
states in GaAs substrate by injecting non-equilibrium 
transient carriers during the pulse duration of 40 fs. Sub- 
sequent carrier-carrier/carrier-phonon collisions lead to 
decoherence and then to a quasi-thermal, hot transient 
carrier distribution within 1 ps. This hot carrier distribu- 
tion further cools down by emitting phonons and eventu- 
ally relaxes back to the equilibrium by carrier recombina- 
tion on the order of 1 ns. In the weak excitation regime 
below 19 /iJ/cm 2 (see Fig. 1), the injected carriers make 
the GaAs regions in the gaps of SRRs and in between 
them conductive. This leads to the reduction of the res- 
onant absorption near uj\ and even complete quenching 
of the resonances [13]. Thus the observed crossover from 
the photoinduced quenching to re-emergence of the THz 
resonance clearly indicate a new collective excitation set- 
ting in at the high excitation regime, not accessible in 
the previous measurements. 

To gain quantitatively insights into the transition den- 
sity N c for the excitation-induced crossover behavior, we 
experimentally determine temporal evolution of the tran- 



sient carrier density. The phase-sensitive nature of the 
time-domain THz field measurement can directly yield 
both real and imaginary parts of the optical conductiv- 
ity in the photoexcited GaAs substrate, on equal foot- 
ing without any model assumptions. This gives a direct 
determination of the transient carrier density based on 
an analysis using the Drude model (see supplementary 
material). Figure 3(a) plots the temporal evolution of 
the transient carrier concentration excited by 70 /iJ/cm 2 , 
which shows an exponential decay with a peak density 
10 17 cm -3 . By correlating the temporal profile to the 
fluence-dependent carrier concentration at Ar =20 ps 
(inset, Fig. 3(a)), one can define a crossover point at 
density N c -4.7 x 10 16 cm -3 at time delay of Ar c -1.73 
ns (black square), corresponding to the transient density 
for 19 /iJ/cm 2 . This naturally divides the transient elec- 
tronic response of photoexcited metamaterials into two 
regimes with opposite time and excitation dependence, 
as illustrated in the shaded areas of Fig. 3(a). Transient 
spectra at various time delays, as shown in Figs. 3(b)-(e), 
show a reversal of the THz transmission change once the 
system is driven cross the boundary, corroborating well 
of our conclusion, e.g., the THz transmission recovers its 
strength at Ar =2.2 ns after initial quenching at time 
delays before 1.73 ns. 

Next, we analyze the transient THz response of pho- 
toexcited metamaterials by full-wave numerical simula- 
tions. Before we take into account any complex processes, 
we first simply consider the photoexcitation-induced dy- 
namic response of the GaAs substrate beneath the SRRs. 
A photoexcited layer is modeled by the well-established 
Drude model, in which, the frequency-dependent com- 
plex conductivity a — Squj 2 / (7 — iuS). Here, So is the per- 
mittivity of vacuum, 7 represents the collision frequency, 

and Up = yj Ne 2 / (eom*) the plasma frequency. Some 
other parameters include the carrier density TV, the free 
electron charge e (1.6 x 10 _19 C) and the effective carrier 
mass in GaAs, m* = 0.067mo (mo, the mass of a free 
electron). Thus, the photoexcited layer has the dielectric 
function e(u) = e s + ia/(eou), where e s (=12.7) is the 
dielectric constant of undoped GaAs. In our simulations, 
we keep 7 = 1.8 THz for all the cases as in Ref. 17. 

Figure 4(a) shows the simulated transmission spec- 
tra for various carrier density N. The simulation clearly 
shows the crossover from a reduction to an enhancement 
of the resonant THz absorption ~1.6 THz as the ex- 
citation from low to high density regime. This agrees 
well with the photoinduced reversible modulation demon- 
strated in the experiments. Specifically, in the low den- 
sity below 4.0 x 10 16 cm -3 , the THz resonance peak in 
the transmission curve exhibit clearly weakened transi- 
tion strength, a blue-shift and spectral broadening. In- 
creasing the carrier density above 4.0 x 10 16 cm -3 , the 
resonant THz peak progressively reemerges, e.g., a clear 
transmission dip seen for the case N = 8.1 x 10 17 cm -3 . 
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FIG. 4. (color online) (a): The simulated THz transmission 
spectra for various carrier densities N. The photoinduced THz 
resonances in the strong excitation regime for two different 
lattice constants 45 and 50 /im: (b) experiment and (c) sim- 
ulation. 



From the simulations, we can attribute the reemergence 
of the resonant THz absorption at high density to the 
effect of first-order diffractive waves from photoinduced 
transient grating. Due to the periodic structure of the 
SRRs, the first-order diffractive wave is found to propa- 
gate in the substrate at high excitation, which strongly 
influence the fundamental absorption mode. This leads 
to a new THz transmission dip at fi = c/(na) above 
a transition density ~ 4.0 x 10 16 cm -3 . Here c is light 
speed in vacuum, n the refractive index of GaAs and a 
the lattice constant of the structure. The SRR structure 
with the dipole resonance close to f\ allows to achieve a 
reversal of resonant THz absorption with increasing pho- 
toexcitation, due to a crossover in the absorption from 
the dipole mode to the diffraction mode. 

To corroborate our model, we compare the re-emerged 
THz resonances from the experiment (Fig. 4(b)) to those 
from the simulation (Fig. 4(c)) for two different lat- 
tice constants, 50 fim and 45 /im, respectively. The 
simulations are performed to the structures with the 
same double- SRR array for the case with the density 
N = 8.1 x 10 17 cm -3 . It is clearly visible that the simu- 
lated transmission dip to higher frequency by decreasing 
the lattice constant, and the dip positions tightly link to 
the frequency of the first-order diffraction mode for each 
case. An extremely good agreement between the simula- 
tion and the experiment is found for the transmission-dip 
position and its shift, which underpins the effect of the 
high-order diffraction mode in the strongly excited meta- 
materials. 

In summary, we have demonstrated ultrafast photoin- 
duced reversible modulation of resonant THz absorption 
in strongly photoexcited metamaterials. Increasing the 



excitation from below to above the threshold density N c , 
we observe a crossover from a complete quenching to a 
reemergence of the THz resonance. Our analysis and 
theoretical simulations, based on the first order propaga- 
tion from the photoinduced transient grating, explain the 
density- and lattice-constant-dependent frequency shift. 
Our results clearly identify the importance of photoin- 
duced optical modes in metamaterials, besides the elec- 
tronic resonances, which should be carefully considered in 
designing future multifunctioning photonic-electronic de- 
vices using the artificial periodical structures. The easy 
implementation of the revealed scheme represents an- 
other approach to achieve nonlinear and frequency- agile 
functions in THz devices. 
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